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ABSTRACT 

In order to find the most extreme dust-hidden high-redshift galaxies, we select 196 extremely red objects in 
the Ks and IRAC bands (KIEROs, 4.5^m]AB > 1-6) in the 0.06 deg- GOODS-N region. This selection 
avoids the Balmer breaks of galactic spectra at z < 4 and picks up red galaxies with strong dust extinction. The 
photometric redshifts of KIEROs are between L5 and 5, with ^ 70% at z ^ 2-4. KIEROs are very massive, 
with M^, ~ 1O'"-1O'^M0. They are optically faint and usually cannot be picked out by the Lyman break 
selection. On the other hand, the KIERO selection includes approximately half of the known millimeter and 
submillimeter galaxies in the GOODS-N. Stacking analyses in the radio, millimeter, and submillimeter all show 
that KIEROs are much more luminous than average 4.5 fim selected galaxies. Interestingly, the stacked fluxes 
for ACS-undetected KIEROs in these wavebands are 2.5-5 times larger than those for ACS-detected KIEROs. 
With the stacked radio fluxes and the local radio-FIR correlation, we derive mean infrared luminosities of 2- 
7 X 1O'^L0 and mean star formation rates of 300-1200 Mq yr~' for KIEROs with redshifts. We do not find 
evidence of a significant subpopulation of passive KIEROs. The large stellar masses and star formation rates 
imply that KIEROs are z > 2 massive galaxies in rapid formation. Our results show that a large sample of dusty 
ultraluminous sources can be selected in this way and that a large fraction of high-redshift star formation is 
hidden by dust. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — radio continuum: galaxies — submillimeter: galaxies 



1. INTRODUCTION 

Deep imaging in the submillimeter opened a new window 
for s t udying high-r edshift ga laxies (ISmail. Ivison. & BlainI 
[19971 iHughes et al.i ri998: Ba rgeret al.lll998h . The brighter 
submillimeter galaxies (SMGs, typical total infrared lumi- 
nosity greater than lO'^L©) that have radio counterparts 
were found to be at redsh ifts mostly between 2 and 3 
jChapman et al.l 120031 l2005h . Despite the extremely large 
star formation rates (SFRs), typically > 1000 Mq yr"', 
such galaxies do not emit strongly in the rest-frame ultra- 
violet (UV) and therefore are gener a lly not Lyman break 
galaxies (LBGs; iPeacock et all 120001: iChapmanetaLl 120001: 



2002 : iKnudsen. van der Werf. & KneibI l2008t IChen et al.1 



Webb et all 



Capak et al, 



20031) . (A counter example is a z = 4.55 SMG in 
12001 which was originally selected as an LBG.) 



This implies a large fraction of dusty star formation at high 
redshift being missed by rest-frame UV surveys. 

This raises two general questions: will we start to see galax- 
ies similar to optically selected ones if we are able to probe 
deeper in the submillimeter, and exactly how much high- 
redshift star formation is missed by optical surveys? We are 
not close to being able to answer either question. Deep sub- 
millimeter surveys in lensing clusters are able to probe dusty 
galaxies with IR luminosities < 10'^Lf:i and constrain their 
number density (iBlain et al.l Il999t ICowie. Barger. & KneibI 
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201 ll) . but the identification of their optical counterparts is 
not yet complete, and the sample sizes are very small. Fur- 
thermore, a z > 4 SMG, GOODS 850-5, was recendy found 
to be extremely faint in the o ptical and near-infrared (NIR) 
at Anh s < 2.2 fim (IWang et al.ll2007l: IWang. Barger. & Cowiel 
120091) . If a significant fraction of the z > 2 SMGs in the IR 
luminosity range lO'^-lO'^ L© are Uke GOODS 850-5, then 
identifying them in the radio and optical will be challenging 
even with next-generation instruments. We therefore seek a 
NIR diagnosis of extremely optically faint SMGs. A gen- 
eral description of the Spitzer Infrared Array Came ra (IRAQ 
colors of bright SMGs was recently developed by lYun et all 
( l2008l) . Here we focus on the reddest galaxies in the IRAC 
and Ks bands. 

In this paper we present a new color selection of extremely 
red objects (EROs) with Ks and IRAC colors of Ks -4-.5 fim 
> 1.6. We refer to such sources as KIEROs. This selection 
is motivated by the fact that at least half of known SMGs are 
redder than this color (Section|2]i. Among all existing ERO se- 
lections, the KIERO selection utilizes the longest wavebands 
that are practically accessible for deep imaging. Unlike some 
other selections for high-redshift red objects, the KIERO se- 
lection does not particularly aim at the 4000 A Balmer breaks 
in galactic spectra, which only enter the Ks band at z ~ 4. 
Instead, this selection aims at galaxies at z > 2 whose ex- 
tremely red colors are likely caused by large dust extinction 
(SectionlHi. Also because of this, we choose the 4.5 /im band 
instead of the 3.6 /im band, to be more sensitive to galax- 
ies whose spectral energy distributions (SEDs) are red over a 
broad wavelength range (cf. sharp spectral breaks). 

The paper is organized as following. We describe the data 
and the KIERO selection in Section |2] the number counts in 
Section|3] the optical and NIR SEDs, redshift distribution, and 
stellar populations in SectionH) and the radio, millimeter, sub- 
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Fig. 1. — Kg, 3.6 fim, and 4.5 /^m color-color diagram. Small gray dots 
are 3 cr 4. 5 fim sources in the catalog of W 10. Spectroscopically identified 
sources in IBarger. Cowie. & Wand j2Q08l) . including stars, are shown with 
color symbols. Radio and submillimeter interferometrically identified 4 cr 
SMGs in the sample of Wang et al. (2004) and Perera et al. (2008) are shown 
with filled squares. The horizontal line indicates our selection criterion for 
KIEROs. 

millimeter, and X-ray properties in Section|5] We identify ac- 
tive galactic nuclei (AGNs) in KIEROs in Section |6] We esti- 
mate the SFRs and SFR densities (SFRDs) of KIEROs in Sec- 
tion |7] We compare KIEROs with other high-redshift galaxy 
populations in Section |8] We discuss the roles of KIEROs 
in galaxy formation in Section |9] and summarize our results 
in Section [TO] Throughout the paper, we assume cosmolog- 
ical parameters of Ho = 71 km s"' Mpc"', = 0.27, and 
f^A = 0.73. All magnitudes are in the AB system, where an 
AB magnitude is defined as AB = 8.9-2.51og(fluxinJy). 

2. THE KIERO SAMPLE 

We use our p ublicly released Ks and IRAC catalogs 
jWang et al.ll2010l hereafter WIO) of the Great Observatories 
Origi ns Deep Surveys-North (GOODS-N; Giavalisco et all 
120041) field. The Ks source catalog of WIO is extracted from 
an extremely deep Canada-France-Hawaii Telescope Ks im- 
age, which has a Icr depth of 0.12 fiSy in the GOODS-N re- 
gion. WIO used the GOODS-N Spitzer IRAC images at 3.6- 
8.0 fim (M. Dickinson et al. 2011, in preparation) to mea- 
sure the IRAC photometry of the Ks selected sources using a 
CLEAN-like method that takes the Ks catalog and image as 
priors. WlO's primary Ks and IRAC catalog contains 15018 
Ks selected sources that are detected at 3.6 and 4.5 fim in 
the GOODS-N region. WlO's secondary catalog contains 358 
/iT^-faint IRAC detected sources. 

We show the Ks, 3.6 fim, and 4.5 /^tm color-color diagram 
for sources with > 3a 4.5 iim fluxes in WlO's primary catalog 
in Figure[T] To obtain colors of GOODS-N SMGs, we looked 
for Ks and IRAC counterparts of the SMGs in the 4 a SCUBA 
sample of Wang. Cowie. & Bargeii ( 120041) and AzTEC sam- 
ple of iPerera et al.. QOQSi). We only include SMGs with un- 
ambig uous radio interferometric identificat ions (Wang et al. 
2004; iPope et al. |l2006>: ' Chapin et al.ll2009j) or submilUmeter 
interferometric id entifica tions ( lono et al. | 120061 : IWang et aLl 
120071: iDaddieTaLiaOOgat I Wang et al.ll201 It A. J. Barger et al. 
2011, in preparation). There are 28 such SMGs and they are 



shown with squares in Figure [T] Fourteen of them are redder 
than Ks -4-.5^m = 1 .6 (horizontal line in the figure) and three 
are not even detected in Ks (lower limits in the figure). They 
occupy a color space significantly different than that of most 
sources detected at both Ks and 4.5 fim. Because mos t spec- 
trosco pically identified SMGs are at z < 3 ( Chapman et al.l 
I2003L E005), the extremely red colors of these SMGs must 
be caused by very high redshifts and/or extremely large dust 
extinctions ( see Se ction 4), like in the case o f GOODS 850-5 
(IWang et al.ll2"007t IWang. Barger. & Cowiell2009h . Such ex- 
tremely red SMGs are of particular interest. Since it is known 
that a variety of optical/NIR criteria are required to include 
all bright SMGs ( iReddv et al.ll2005h . here we only focus on 
selecting and understanding the extremely red ones. 

In order to find sources similar to the extremely red SMGs 
but slightly fainter in the submillimeter and therefore unde- 
tected by current submillimeter surveys, we select KIEROs 
with Ks -4-.5 fim > 1.6. We require that the sources be de- 
tected at > 3a in the 4.5 band and in at least one other 
IRAC band. For galaxies detected in the IRAC bands but not 
in the Ks band (in WlO's secondary catalog), we adopt their 
la Ks limits. Although WlO's primary catalog is a Ks selected 
one, the above inclusion of /Js-faint 4.5 /.tm sources in WlO's 
secondary catalog makes the selection of KIEROs in principle 
a 4.5 fim selection. We visually inspected the IRAC images of 
all such galaxies and excluded those whose IRAC fluxes are 
less reliable, mainly galaxies affected by very bright stars in 
the field and galaxies blended with > 3 bright nearby galax- 
ies. These are the limitations of WlO's catalogs (see WIO for 
discussion). This reduces the number of selected KIEROs by 
< 10% and should not impact our analyses, even if some of 
them are truly red sources. We are left with 104 KIEROs in 
the main /T^-selected sample and 92 in the Ks undetected sam- 
ple. Among the above mentioned 28 identified SMGs in the 
GOODS-N, 14 are KIEROs. 

3. NUMBER COUNTS 

Figure |2] shows a Ks and 4.5 /im color-magnitude diagram 
for KIEROs and field galaxies. The KIEROs are a relatively 
rare population among all sources detected at 4.5 /im, and 
the selection of KIEROs lies very close to the detection limits 
at both Ks and 4.5 ijm. At Ks, WIO provided a description 
of the completeness limits (e.g., the Ks 30% limit shown by 
the diagonal dashed line in Figure |2|i. On the other hand, it 
is nontrivial to quantify the completeness at 4.5 iim, which 
depends on both the Ks completeness and the very complex 
procedures used by WIO to construct the 4.5 /im source lists. 

Because of the above issues, we have decided to present 
the Ks number counts of KIEROs (Figure|3]l, even though this 
is a 4.5 /im selected population. In our survey there are no 
KIEROs brighter than 6 fiiy at Ks- The raw counts slightly 
flatten at ^ 1 /iJy and drop rapidly below 0.3 /iJy, consis- 
tent with the Ks detection completeness. Without consider- 
ing the completeness at 4.5 /im, simply correcting the counts 
with the Ks completeness derived by WIO suggests that the 
slope of the counts does not change significantly down to a 
flux level of - 0.3 /iJy. A slope of 0.41 ± 0.07 in the logN- 
mag space was fitted at Fks > 0.3 /iJy using completeness 
corrected counts weighted by the Poisson errors. Comparing 
to the slope for the entire Ks sample in the same flu x range 
(- 0.2, e.g.. lMaihara et alllioOll: iKeenan etal]|2010h . this is 
much steeper. This is likely a consequence of the higher red- 
shifts of KIEROs (see next section). Furthermore, the large 
number of AT^-undetected sources that meet our KIERO selec- 
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Fig. 2. — Ks and 4.5 color-magnitude diagram. Dots are field galaxies 
detected (3(t) in both the Kg and 4.5 fj,m bands, as well as at least one other 
IRAC band. Diamonds are Ks detected KIEROs, and triangles are Icr color 
lower limits of Ks undetected KIEROs. Squares are the radio and submil- 
limeter interfeormetrically identified SMGs. The horizontal solid line shows 
the selection criterion for KIEROs. Some of the galaxies above the horizon- 
tal line are not designated as KIEROs because their 4.5 fim fluxes in WlO's 
catalogs are less reliable, mostly bad photometry affected by very bright stars 
in the field or galaxies blended with > 3 bright nearby galaxies. The diagonal 
dashed line shows the Ks 30% completeness limit determined by WIO. The 
vertical dotted line shows the 3a limit of the 4.5 /-tm fluxes from WIO. The 
completeness at 4.5 fim is unclear due to the nature of the 4.5 fim catalog 
(see WIO for details). 
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Fig. 3. — Ks number counts of KIEROs. The solid histogram shows the 
number of sources detected in the Ks band at > 3a in bins of 0.2 dex in 
flux. The dashed histogram shows the same for sources with only Ks upper 
hmits (Icr). The vertical dotted lines are the 90% and 30% completeness 
limits of the WIRCam Ks image derived by WIO. The actual completeness 
of the Ks- 4.5 /im selection should be even lower (see text). The number 
distributions shown here are not corrected for completeness. 

tion criterion (dashed histogram in FigureO suggests that the 
counts may be still rising at a Ks flux of ^ 0.2 /iJy. 

4. REDSHIFT AND STELLAR POPULATION 

To obtain a rough idea about the SEDs and redshifts of 
KIEROs, we plot the color-redshift diagram of spectroscop- 
ically identified 4.5 /im sources in Figure |4] Most 4.5 /im 
sources have 7*^5-4.5 /im colors bluer than 1.0. Even AGNs, 
which tend to have unusual colors, are mostly bluer than 1 .0. 
In general, we only expect heavily extinguished (Ay > 2) 
galaxies at z > 2 to be redder than 1.6. At z < 2, galaxies 
can be redder than 1.6 in only the most extreme cases. We do 
not expect to see many such galaxies. To verify that KIEROs 
are mostly at z > 2, we looked at the spectroscopic redshifts 
of GOODS-N galaxies and carried out a photometric redshift 
analysis. 
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Fig. 4. — Ks-4.5 fj,m color vs. spectroscopic redshifts. Dots are 3a 4.5 ^m 
sources in the catalog of WIO that are spectroscopically identified in Barger 
et al. (2008). Squares show X-ray selected AGNs (soft or hai'd X-ray lumi- 
nosities greater than lO**^ erg s"'). Curves show colors derived from the mod- 
els of elliptical galaxies (thin red curves ICol eman, Wu, & Weedman 1 98(3) 
and starburst galaxies (thick blue curves, iKmnev et al.iil996i) . reddened with 
the Calzetti et al. (2000) extinction law with A v = (solid), 1.0 (dotted), 2.0 
(dashed), and 3.0 (dash-dotted). 

4.1. Spectroscopic Redshift 

In the entire KIERO sample, only two sources are bright 
enough in the optical a nd NIR to have spectroscopic re dshifts 
in our redshift survey dBarger. Cowie. & Wang||2008l) . Both 
are detected in the 2 Ms C handra image (source numbers 109 
and 135 in the catalog of [Alexander et al.]l2003h . One is the 
SMG GOODS 850-7 (F1.4 ghz = 58 fiJy, Fgso = 6.2 mJy) 
at z = 2.578. The other is a F14 ghz = 83 ^Jy radio source at 
z = 1 .790. This radio source is very close to the bright SMG 
GOODS 850-2 {F^5u = 10.3 mJy), with an angular separa- 
tion of ~ 8 "6 based on its SCUB A position (Wang et al. 2004; 
see also lBarger. Cowie. & Richards 2000) . However, our re- 
cent submillimeter interferometric observation (A. J. Barger 
et al., in preparation) shows that GOODS 850-2 is another 
faint Ks-hand source that also enters our KIERO sample. 

In addition to the two optical spectroscopic redshifts, an 
SMG KIERO, GOODS 850-5 (Wang et al...2007) . has a plau- 
sible millimeter spectroscopic redshift of 4.042 based on a 
line detection at 9 1 .4 GHz that is likely a CO(4-3) transition 
dDaddi et al.ll2009bh . The three spectroscopic redshifts span a 
redshift range that agrees with what would be expected based 
on Figure |4] However, the sample size is too small and we 
need to rely on photometric redshift fitting. 

4.2. Photometric Redshift 

We carried out a photometric redshift analysis on a subsam- 
ple of 76 high S/N KIEROs whose Ks fluxes are greater than 
0.2 /iJy. We included the photometry in the U, HST ACS, 
J, WFC3 F140 W, Ks, and IRAC ba nds. We adopted the U- 
band fluxes from lCapak et al.l (120041) . For the ACS bands, we 
directly ado pted the fluxes from the GOODS-N v2 catalog 
( iGiavalisco et al. 20041). We measured 7-band fluxes from a 
CFHT WIRCam image that contains ^ 27 hr of total integra- 
tion. We obtained the /-band images from the public archive. 
They were originally obtained by a Taiwanese team led by 
Lihwai Lin (2010, in preparation). We reduced and processed 
the y-band images in an identical manner to WIO for the Ks- 
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-0.01 in Az. We also made the same comparison on Hyperz 
results. The outlier fraction and dispersion in Az are both 
significantly worse than the EAZY ones. 
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Fig. 5. — Results of photometric redshifts for Fa:., > 0.2 fi]y sources in the 
KIERO sample. A total of 76 KIEROs are included. The top panel shows 
the distribution of the EAZY photometric redshifts. The bottom panel shows 
the stellar masses fi tted by Hyperz at t he EA ZY redshifts using the models of 
IBruzual & Charloll t2003l) (see Section l4!3t . The vertical bars show the max- 
imum and minimum stellar masses from various types of models, while the 
squares are the best-fit masses. Large open squares denote AGNs identified 
in Section|6] which may have problematic photometric redshifts and masses. 



band images. The HST WFC3 F140W data will be described 
by A. Barger et al. (2011, in prep). 

For photometric redshift fitt ing, we tried two packages; 
the la test version of Hyperz dBolzonella. M iralles. & Pel^ 
|2000B, and EAZY dBrammer. van Dokkum. & Co ppi 2008). 
We found that EAZY generally produces better results and in- 
cludes a necessary feature for our studies (see below). Our 
primary photometric redshift results in this paper are thus 
b ased on EAZY. We adopted the default set of SED templates 
of lBlanton & Rowelsl ( l2007h . provided in the EAZY package. 
This set includes five templates ranging from very blue and 
young galaxies with strong emission lines, to galaxies dom- 
inated by old stars. We included an additional SMG-type 
dusty starburst model (f = 50 Myr, Ay = 2.75), also provided 
in the EAZY package. We refer to Brammer et al. (2008) 
for detailed descriptions on all these SED models. The tem- 
plates here all already include certain amounts of extinction 
in order to fit galaxy SEDs in deep surveys. To account for 
very dusty sources in the KIERO sample, we further red- 
dened these templates (including the already reddened SMG 
template) by Ay = , 0.5, and 1.0 with the extinction law of 
ICalzetti et al.l (120001) . To account for the situation where there 
are more than two distinct stellar populations in a galaxy, we 
allowed for all linear combinations of the templates in the fit- 
ting. In Section |43] we will show that the EAZY feature of 
allowing for all linear combinations of templates is essential. 
Because a fitted SED is a combination of multiple templates, 
we are not able to quote an extinction value for a source. In 
Figure |3a) we show the distribution of the EAZY photomet- 
ric redshifts. 

In order to see how reliable the fitting is, we ran EAZY on 
^ 1500 Fks > 0.2 /xJy galaxies in the spectroscopic sample 
of Barger et al. (2008). The results are shown in Figure |6l 
If we define Az as (Zph-z.sp)/(l +z.sp), then only 5.3% of the 
galaxies have | Az| > 0.2. After excluding these outhers, there 
is an rms dispersion of 0.05 and a median systematic offset of 

* see also |http://www.ast.obs-mip.fr/users/roser/hyperz7] 
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Fig. 6. — Results of EAZY photometric redshifts for the Fks > 0.2 /iJy 
sources (small dots) in the spectroscopic sample of Barger et al. (2008). A 
total of 1987 galaxies are included. Large squai'es are the two KIEROs with 
spectroscopic redshifts. 

Although the results in Figure|6]look good, there is no guar- 
antee that we can achieve the same on KIEROs. The two 
KIEROs with spectroscopic redshifts at z = 2.578 and 1.790 
have EAZY photometric redshifts of 2.03 and 1.96, respec- 
tively (large squares in Figure|6l). These are excellent, but the 
sample is too small for us to comment on the overall quality. 
We can also compare the photometric redshifts on KIEROs 
that show continuum dropouts (Lyman breaks, see also Sec- 
tion l8.3l ). There are 15 fe-dropouts in the KIERO sample, 6 of 
which are bright enough at Ks to have photometric redshifts. 
Their mean photometric redshift is 3.72 ± 0.50. There are 4 
v-dropouts in the KIERO sample, 2 of which have photomet- 
ric redshifts of 5.02 and 4.71. These are all consistent with 
the ex pected redshifts for the dropouts (e.g., iBouwens et alj 
l2007h . Finally, Figure |7] shows that the majority of the fit- 
ted SEDs reasonably represent the observed SEDs, including 
various spectral breaks (the Lyman breaks, Balmer breaks, 
and 1 .6 fim bumps), except AGNs (Section |6]l with feature- 
less power-law SEDs in the IRAC bands. It is worth noting 
that there are three non-AGN data points at z > 4 with un- 
usually large stellar masses (Section l43] ) of ^ lO'^M©. They 
are #23, 24, and 38 in Figure |7] Their observed SEDs are 
also quite featureless in the IRAC bands, and all fade dramat- 
ically at Kg and shorter wavebands, likely because of extinc- 
tion. Although they are not identified as AGNs, we suspect 
that their photometric fittings are unreliable. Except for AGNs 
and these few problematic cases, the EAZY photometric red- 
shifts for t he m ajority of the sources seem reasonably good. 
In Section 15.21 we present independent evidence suggesting 
that the redshift distribution in Figure [Sja) is approximately 
correct. 

We did not attempt to study the redshifts of fainter sources 
(Fks < 0.2 /iJy). Such sources often have robust photometry 
in only two IRAC bands, which is insufficient for good photo- 
metric redshifts. One might expect that these fainter sources 
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Fig. 7. — Observed (diamonds) and fitted (curves) SEDs in the EAZY pliotometric redsliift analysis. Sources witliout symbols in the ACS and F140W bands 
ai'e outside the image region. Sources without symbols in the J band do not have reliable /-band counterpart identifications. AGNs identified in Section|6]are 
indicated in the upper-left corners and have problematic photometric redshifts. 



should have higher redshifts rather than lower intrinsic lumi- 
nosities. However, we cannot test this with the current data. 

4.3. Stellar Population 

When testing the photometric redshift fitting, we found that 
many KIEROs have SEDs that cannot be fitted with a sin- 
gle stellar population by Hyperz. For example, on the above 
mentioned six /7-dropouts, Hyperz returns a mean redshift of 
2.98 ±0.34. This is clearly too low, a consequence of trying 
to fit the complex SEDs with single stellar populations. This 
is a key reason for our adoption of EAZY, which allows for 
combinations of different SED templates. The most obvious 
cases of such KIEROs are sources 3, 11, 30, 35, 41, 53, and 
70 in Figure |7] In such KIEROs, there are well-evolved old 
stellar populations, which produce the observed NIR lumi- 
nosity, as well as unobscured ongoing star formation, which 
produces the strong rest-frame UV emission. This property 



is similar to that of o bjects selected with z-3.6 i-im > 3.25 
(lEROs) bvlYan et al.l(l200l in the Hubble Ultra Deep Field 
(HUDF: iBeckwith et al.l 20061) . Interestingly, our stacking 
analyses in the radio and FIR (Section |5] and Table [T]) show 
that sources that are faint in the rest-frame UV have more 
intensive starbursts. SFR estimates based on rest-frame UV 
may miss the strong dust-hidden star formation on such opti- 
cally faint sources. 

Another interesting question to ask is whether or not 
KIEROs are massive galaxies. Unfortunately, EAZY does not 
provide direct estimates of stellar masses. To do this, we rely 
on Hyper z, which builds in the ste llar population synthesis 
models of lBruzual & Charloll (l2003h and provides stellar mass 
estimates. We forced Hyperz to fit at the EAZY redshifts. We 
derived two sets of stellar masses by fitting to the full SEDs 
of the KIEROs, and by fitting to just the infrared SEDs in the 
J and redder bands. We found that neither of the two masses 
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Fig. 7. — Continued. 



systematically favor high or low mass values. For roughly 2/3 
of the objects, these two masses agree within 0.4 dex. We thus 
believe that for most of the objects, the derived masses are ro- 
bust. On the other hand, since many of the KIEROs show 
two different stellar populations, it is better to just fit the rest- 
frame NIR parts of the SEDs. The mass-to-light ratios in the 
NIR are relatively insensitive to both extinction and ongoing 
star formation. Nevertheless, we remind the reader that our 
mass estimates are significantly limited by the fact that Hy- 
perz can only fit the SEDs with single stellar populations. 

In Figure |5lb) we show the best fit stellar masses (squares) 
and the maximum and minimum masses (v e rtical bars) from 
the various templates of iBruzual & CharloB ( l2003h . based on 
just fitting the infrared SEDs. For most KIEROs, the data only 
allow small ranges of stellar masses, a consequence of similar 
mass-to-light ratios in the NIR. They have large stellar masses 
of 10'"Mq to > 1O'^M0. Some of the masses seem unusually 
large, especially at the high-redshift end {z > 4) where less 



photometric data points are available and there is a degeneracy 
between photometric redshift, age, and extinction. This is a 
fundamental limit of the current data. We also note that there 
are 18 sources with stellar masses less than IO^^Mq, giving 
a mass distribution that is somewhat disjoint from that of the 
other KIEROs. These are all optically bright sources with 
^F850LP > 0.2 i^Jy and have lower redshifts. They are likely a 
less dusty subclass of KIEROs. 

In summary, the majority of KIEROs in the photometric 
redshift subsample are more massive than 1O'°M0. In ad- 
dition to the massive stellar populations, a large fraction of 
KIEROs show young stellar populations in the rest-frame UV. 
In the following sections of this paper, we will focus on the 
star formation properties of KIEROs. 

5. MULTIWAVELENGTH PROPERTIES AND STACKING ANALYSES 
5.1. Radio Properties 
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In order to understand the star formation activities in 
KIEROs (albeit with AGN contamination), we studied their 
radio properties. We used the Very Large Array (VLA) 
1 .4 GHz image and ca talog in the Hubble D eep Field-North 
(HDF-N) pubHshed bv lMorrison et all (120101) . The image has 
a 5 (T sensitivity of 20 /iJy at the field center, and the catalog 
contains ^ 430 sources in the GOODS-N regi on. Among the 
196 K IEROs, 22 are included in the catalog of lMorrison et alj 
(l2OT0h . 

For sources not in the catalog of lMorrison et al.l ( 120101) . we 
measured their fluxes from the image with 3" apertures af- 
ter correcting for the VLA primary beam and the band-width 
smearing. We compare d our fluxes on c ompa ct sources with 
those in the catalog of iMorrison et al.l (120101) . and they are 
fully consistent with each other The adopted aperture size 
is 2 X the VLA beam FWHM and corresponds to 20-25 kpc 
at z = 2-5. This size is slightly less optimized for detecting 
faint and compact sources, but it encloses most of the radio 
flux if the radio-emitting region in a high-redshift galaxy is 
offset from the position determined by NIR observations by 
a few kpc. Such offsets are not uncommon in nearby inter- 
acting/merging systems. The measured 1.4 GHz fluxes of the 
KIEROs are shown in the bottom and middle panels of Fig- 
ure H] 

AB Magnitude 
25 24 23 22 21 20 
^ .jQQQ^- ■ 53j±2.1 1 

>- 21.5±1.5 ' ; 



3 10.0- 6.8±0.9 " ^ 




1. 400 : 



I 200 - H - 



1 10 

F4.5 (iiJy) 

Fig. 8. — Radio flux vs. 4.5 /xm flux for KIEROs. The bottom and mid- 
dle panels show the radio fluxes of individual sources. The middle panel is 
a blowup of the bottom panel. Sources in the bottom and middle panels en- 
closed in large squai'es are AGNs identified in Section|6] The top panel shows 
the stacked radio fluxes in each 0.5 dex 4.5 /^m flux bin. Solid squares are the 
mean fluxes after excluding the Fi 4 ghz > 600 /iJy bright source. The exact 
values are also shown near the symbols. Open squares are the median fluxes, 
which further avoid influence from > 200 /^Jy sources. Triangles show the 
mean fluxes of the entire 4.5 ^m sample in the same flux ranges. The dotted 
line corresponds to roughly 5a. Stacking of faint sources below this line is 
presented in the lower half of Table[T] 

We can study the mean 1 .4 GHz fluxes of KIEROs via a 
"stacking analysis," in which we measure radio fluxes with 3" 
apertures at the locations of KIEROs and average the results. 
First, we need to understand the bias and uncertainty in such 
stacked fluxes, for which we took a Monte Carlo approach. 



We randomly placed 3" apertures in the VLA image, calcu- 
lated their mean flux, and repeated a large number of such 
measurements. The mean flux measured in this way is an 
estimate of contamination from physically unrelated nearby 
sources. We refer to this as confusion. We excluded flux val- 
ues greater than 600 to avoid being biased by the very 
bright source, which is likely an AGN (also see Section l53i l. 
In the IRAC region, the mean 1 .4 GHz flux is 0.41 /iJy for the 
random positions in the Monte Carlo simulations. We sub- 
tracted this value from the mean flux of each KIERO sample. 

We treated the dispersion between various measurements 
of the mean flux of random apertures as noise in the mean, 
which has contributions from the receiver noise and from the 
confusion noise. We varied the number of random apertures 
that we placed in the image each time, and we found that even 
with source densities up to 10^ deg~^, the measured noise in 
the mean precisely scales with the square-root of the source 
number In the IRAC region, such noise is 8.7 /iJy per source. 
This value is larger than the nominal 4 /^Jy noise of the image, 
a consequence of confusion caused by bright sources. On the 
other hand, the average fluxes we determine are considerably 
below the conventional single-source detection limit. 

We present our major stacking analysis results in Table [T] 
For the entire sample of 196 KIEROs, excluding the > 600 
/^tJy one, the mean 1.4 GHz flux is 9.9 ±0.6 ^liy, and the me- 
dian is 3.07 /xJy. For comparison, in the 4.5 /im sample there 
are 8977 sources with 5 > 1 /^Jy- Their mean 1 .4 GHz 
flux is 4.7 ±0.1 /iJy, and their median flux is 2.46 /iJy. Af- 
ter considering the fact that most of the 4.5 sources are at 
z<2 while most of the KIEROs are at z > 2, it is clear that 
our ^^5-4.5 /xm > 1.6 selection picks up a much more radio 
luminous population. 

The high S/N of the mean 1.4 GHz flux for KIEROs al- 
lows us to further break down the sample. In the upper panel 
of Figure [8] we show the mean 1.4 GHz flux of the KIEROs 
as a function of 4.5 /im flux (filled squares). (For compari- 
son, we show the same for the 8977 5 ^,,1 > 1 A^Jy sources 
with open triangles.) There is a strong linear correlation be- 
tween the mean radio flux and the 4.5 /im flux. In our KIERO 
sample the mean 3.6 /im-4.5 /im color is 0.19, which cor- 
responds to a spectral index of a = 0.78 (f^, oc v'"). This is 
similar to the radio spectral index of normal galaxies (0.7- 
0.8), implying similar /iT-corrections in the radio and 4.5 /im. 
Therefore, the correlation could suggest an approximately lin- 
ear relationship between average SFR measured from the ra- 
dio power and stellar mass measured from the rest-frame NIR. 
However, in Section |2l we show that the mean radio fluxes in 
two redshift bins between z = 2 and 4 do not obviously depend 
on the stellar masses. So the linear relationship between radio 
and 4.5 /im fluxes observed here is more likely to be just an 
apparent effect of redshift. 

In Figure |9] we show the stacked radio images of various 
KIERO subsamples. We start by discussing the top row, 
where we have applied a limit of 4 ghz < 600 /^Jy. The 
first thumbnail shows the stack of all the Fi 4 qhz < 600 /iJy 
sources. We then stacked the radio fluxes of the photomet- 
ric redshift subsample of 76 KIEROs. These sources have, 
on average, higher 1 .4 GHz fluxes than other KIEROs, since 
they are selected to be brighter in the NIR, and the radio flux 
correlates with the NIR flux. We find that the KIEROs at 
z < 3 (second thumbnail) and the KIEROs at z > 3 (third 
thumbnail) have similar mean radio fluxes. The observed 
radio flux for a given source has a redshift dependence of 
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TABLE 1 

Stacking Results in the Radio. Millimeter, and Submillimeter 



Samples 


1.4 GHz Stacking 


(0.06 deg2) 


1 100 fim Stacking 


(0.06 deg2) 


850 Stacking (0.028 deg-) 








EBL" 




(F) 


EBL" 




(F) 


EBL" 






(MJy) 


(mJy deg"^ ) 




(raJy) 


(Jy deg"-) 




(mJy) 


(Jy deg"-) 


^4.5 Mm > 1 A^Jy 


8977 


4.7 ±0.1 


706 ±15 


8984 


0.090 ±0.01 


13.4±1.8 


4551 


0.14±0.04 


22.7 ±6.4 




195 


9.9 ±0.6 


32.2 ±2.0 


196 


0.50 ± 0.08 


1.64 ±0.26 


87 


1.44 ±0.28 


4.47 ±0.88 


KIERO, z < 3.0 


38 


22.7 ±1.4 


14.4 ±0.9 


38 


0.56 ±0.18 


0.36±0.12 


18 


2.48 ±0.62 


1.59 ±0.40 


KIERO, z > 3.0 


37 


21.0±1.4 


12.9 ±0.9 


38 


1.57±0.18 


1.00 ±0.12 


14 


3.19±0.71 


1.59±0.35 


KIERO, ACS 


112 


6.0± 1.0 


11.2± 1.5 


112 


0.23±0.11 


0.42 ±0.20 


52 


0.71 ±0.37 


1.32 ±0.68 


KIERO, non-ACS 


61 


16.2±1.1 


16.7±1.1 


62 


0.94±0.14 


0.97±0.14 


32 


2.59 ±0.47 


2.96 ±0.54 




Fi 


4 GHz < 20 fily Sources 


Fi 


100 fim < 3 mJy Sources 




^850 ,jm < 6 mJy Sources 


^4.5 urn > 1 A^Jy 


8282 


1.76 ±0.07 


242 ± 10 


8849 


0.062 ±0.011 


9.2± 1.6 


4446 


0.092 ± 0.037 


14.5 ±5.8 


KIERO 


165 


2.6 ±0.5 


7.14± 1.45 


177 


0.14 ±0.08 


0.41 ±0.23 


78 


0.65 ±0.28 


1.82 ±0.77 


KIERO, z < 3.0 


27 


2.4± 1.3 


1.07 ±0.59 


34 


0.20±0.18 


0.12±0.10 


15 


0.90 ±0.63 


0.48 ±0.34 


KIERO, z > 3.0 


23 


8.6 ±1.4 


3.31 ±0.54 


28 


0.72 ±0.20 


0.33 ±0.09 


9 


0.78±0.81 


0.25 ±0.26 


KIERO, ACS 


102 


1.1±0.7 


1.93± 1.14 


107 


0.09 ±0.10 


0.16±0.18 


49 


0.31 ±0.35 


0.54 ±0.61 


KIERO, non-ACS 


48 


5.0±1.0 


4.00 ±0.80 


51 


0.21 ±0.14 


0.18±0.11 


26 


1.16±0.48 


1.08 ±0.45 



" W is the number of sources included in each sample at the corresponding wavelength. 
" EBL is the surface brightness contribution to the extragalactic background hght. 

In the radio stacking, we always exclude sources brighter than 600 /iJy to avoid bias from radio AGNs (see text). Such bright sources bias the 
stacked fluxes of both the 4.5 ^m population the KIERO population by ~ 30%. 
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Fig. 9. — Mean radio images of the KIERO subsamples in Table[T] The 
top row shows the results of F[ 4 qHz < 600 fi]y sources, and the bottom row 
shows the results of Fi 4 qhz < 20 fi]y sources. From left to right are all 
sources fainter than the above flux limits, z < 3 sources, ; > 3 sources, ACS 
detected sources, and ACS undetected sources. Each panel has a size of 20". 
The small circles indicate our 3" flux apertures. The grayscale ranges in all 
panels are -1 .5 to 3 fiJy beam"' . 

■Pradio {I + zY'" / diizf' , where di is the luminosity distance. 
Using this relation and the measured mean radio fluxes, we 
find that KIEROs at z > 3 are 2 x more radio luminous than 
KIEROs at z < 3 (see also, Section|7]and Table|2]i. This may 
be because the 4.5 fim selection is only sensitive to more lu- 
minous systems at higher redshifts. 

We also stacked the radio fluxes of the ACS-detected 
(fourth thumbnail) and the ACS-undetected (fifth thumbnail) 
KIEROs. There is a striking anti-correlation between the 
radio fluxes and the optical fluxes of KIEROs. The ACS- 
undetected subsample is, on average, 2.5 x brighter than 
the ACS-detected subsample in the radio. The same anti- 
correlation is also clearly observed when using 850 and 
1100 fim fluxes rather than radio fluxes (see next section). 
This shows that the most active star forming galaxies in the 
universe are deeply hidden in dust and can hardly be traced 
with rest-frame UV observations, even with the extremely 
deep ACS observations in the GOODS-N. 

It is also important to ask whether the above results are 
general properties of the entire KIERO population, or are bi- 
ased by a few bright sources. In order to investigate this, 
we removed sources brighter than 20 ijJy in the radio (dot- 
ted line in Figure [8]i and repeated the same stacking proce- 
dures. We adopted the same 20 fiiy cut in the Monte Carlo 
determinations of the confusion effect. The stacking results 
are listed in the lower half of Table [T] As expected, all the 
mean fluxes decrease, but they are still greater than at least 



1.5(7. The general pattern is similar: the KIERO population 
is, on average, brighter in the radio than the 4.5 fim sources; 
and the ACS-undetected KIEROs are brighter in the radio 
than the ACS-detected KIEROs. We show some of the re- 
sults in the bottom row of Figure |9] where we have applied 
the /^L4 GHz < 20 fiJy limit. The thumbnail images corre- 
spond to all of the KIEROs below this limit (first), only the 
z < 3 sources (second), only the z > 3 sources (third), only the 
ACS-detected sources (fourth), and only the ACS-undetected 
sources (fifth). 

In Figure|8]we see that nearly all (85 out of 86) sources with 
^4.5 fim < 1 /^Jy are below the 20 fiiy radio flux cut. Their 
stacked radio flux is 0.42 ±0.76 ^Jy, or 0.18 ±0.73 ^Jy after 
removing the one > 20 /iJy source. Both are consistent with 
a null detection. On the other hand, if we look at sources with 
^4.5 > 1 fJ'iy, the positive signal in their radio fluxes is very 
clear In the 5 = 1-3.2 fi3y bin, the mean radio flux is 
6.76±0.90 jiJy for 65 sources, or3.44±0.88 after excluding 
the six > 20 fiJy sources. 

We conclude that the majority of the NIR bright (F4 5 > 
1 /iJy) KIEROs are luminous in the radio, especially the op- 
tically faint sub class. This is true even if we remove radio 
detected KIEROs from the sample. On the other hand, the 
NIR faint {F4.5 < 1 A*Jy) KIEROs (43% of the total KIERO 
population) do not have detectable radio emission. In Sec- 
tion |7] we will use the stacked radio fluxes to study the star 
formation properties of KIEROs. 

5.2. Millimeter and Submillimeter Properties 

There exist numerous millimeter and submillimeter con- 
tinuu m surveys and follow-up studies in the G OODS-N 



(e.g.. iHughes et all [T99I iBarger. Cowie. & Richards. ,2000t 
Borvs et al.ll2003t Wang et al. 2004; IWang Cowie. & Barged 



2006t iPope et al. 1 120061; Wera et al.l'12008^ iGreve et al.ll2008t 
Chapinet al.1 120091) . Here we considered two source cata- 
logs, our SCUBA 850 /im catalog (Wang et al. 2004) and 
the AzTEC 1100 ^m catalog (Perera et al. 2008) both with 
nearly complete multiwa velength identifications (iPope et al.l 
I2OO6I; IChapin et aLll2009l) . To be conservative, we only con- 
sidered robust sources (4a) with unambiguous 1 .4 GHz iden- 
tifications that we agree with, or with submillimeter interfer- 
ometric identifications. 
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There are 28 850 /im and 1100 /im selected sources with 
such identifications in the IRAC area, 14 of which pass our 
Ks-A-.5 fim > 1.6 criterion. In addition, we found that sev- 
eral unidentified millimeter and submillimeter sources have 
nearby KIEROs within a few arcsec, sometimes multiple 
ones. The high fraction of submillimeter sources is consis- 
tent with our 1 .4 GHz finding that we are selecting active star 
forming galaxies at high redshift with our KIERO selection. 

We also stacked the millimeter and submillimeter fluxes of 
KIEROs. Here we used the 850 /im image of Wang et al. 
(2004) and the 1100 /im image of Perera et al. (2008). The 
stacking is similar to that in the radio. However, the caveat 
here is a possible upward bias caused by the coarse beams in 
the millimeter and submillimeter If the stacked population 
has strong clustering at scales smaller or comparable to the 
beam, then the stacked flux would be p ositively biased. This 
was discussed by Wang et al. (2006) and Serieant et al. (20()8|) 
and recently demonstrated by Marsden et al., (2009). Mars- 
den et al. performed a stacking analysi s using the deep Very 
Large Telescope (VLT) NIR sample of iGrazian et al.l (l2006al 
5 X 10^ sources per deg-) in the GOODS-S and the BLAST 
250-500 /im images (beam size 0'.6-l'). Their stacked far-IR 
(FIR) fluxes from the VLT NIR galaxies exceeded the FIR ex- 
tragalactic background light (EBL) by factors of 2.5-30. The 
authors attributed this to small-scale clustering of the NIR se- 
lected galaxies, which shows a ^ 100% excess in the number 
of galaxies at scales of ^ 1'. 

We performed a similar clustering analysis and found that 
our KIERO population only shows a 7% and 9% excess of 
sources within 15" (SCUBA beam FWHM) and 18" (AzTEC 
beam FWHM), respectively. This is much lower than the 
^ 100% excess of the VLT NIR sample at the size scales of 
the BLAST beams. Thus, we conclude that stacking analy- 
ses of the GOODS-N SCUBA and AzTEC images using the 
KIERO population should not be biased by more than a few 
percent by small-scale clustering. (For comparison, we made 
similar tests on the entire Ks-hand sample. There is a ~ 30% 
and ~ 40% excess of sources at size scales of 15" and 18", 
respectively.) After integrating the excess amounts convolved 
with the beam profiles, we found that these would bias the 
stacking results by up to 10%. This is likely an upper limit, 
since not all clustered objec ts are millimeter or su bmillimeter 
sources, as pointed out by iSerieant et al.l (120081) . We there- 
fore conclude that bias caused by small-scale clustering in our 
sample is unlikely to affect our millimeter and submillimeter 
stacking results. 

Instead of aperture photometry, the 850 and 1 100 fim fluxes 
were measured with weighted point-source filters. The filter- 
ing and weighting schem es are described in Wang et al. (2004) 
and lPerera et al.l i l2008l) . We measured the noise (sky, instru- 
mental, and confusion noises) in the weighted fluxes and the 
flux bias (caused by random nearby sources) in a way iden- 
tical to that for our radio stacking. Our noise measurements 
are 2.99 and 1.13 mJy per source for SCUBA and AzTEC, 
respectively. The flux biases are both zero, because of the 
zero sums of the beams. In our SCUBA image, there are two 
negative 50% sidelobes produced by the secondary chopping, 
producing a beam that has a zero total power The PSF of 

* Our group (Wang et al. 2006) did not find a stacked SCUBA 850 ^tm 
flux that exceeds the EBL Hmit, although we used extremely deep and dense 
optical and NIR samples. The difference between our results and that in 
[Marsden et al. (2009) could be attributed to the dramatic difference in beam 
sizes of SCUBA and BLAST, or issues in the methodology adopted by either 
group. 



AzTEC is more co mplex but the filtering function adopted by 
iPerera et"an (l2008h also has a zero sum. 
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Fig. 10. — Millimeter and submillimeter fluxes vs. 4.5 /-tm flux for 
KIEROs. From the top we show stacked 1 100 fim fluxes, 1 100 fim fluxes of 
individual sources, stacked 850 fim fluxes, and 850 fim fluxes of individual 
sources. Individual sources enclosed in large squares are AGNs identified in 
Section|6] In the panels for stacked fluxes, solid squares show error-weighted 
stacked fluxes within bins of 1 dex in 5 for the KIEROs. The exact 
values are shown neai' the symbols. Open squares show median fluxes for the 
same KIEROs in the same bins. For comparison, we show with downward 
triangles stacked fluxes from the 4.5 fim selected sample in 0.5 dex bins. Un- 
like the case in Figure [H these are upper limits and may be overestimated 
by a few to ~ 10 percent because of small-scale clustering (see text). The 
4.5 fim population is not detected at both 850 and 1100 fim at F^ j j,,^ < 1 
/iJy, despite the great numbers of sources there. The dotted lines con'espond 
to roughly Scr. Stacking results for faint sources below these lines are pre- 
sented in the lower half of Table[T] 

In Figure [To] we present the measured millimeter and sub- 
millimeter fluxes of the individual sources in the second and 
fourth panels, respectively, and the stacked fluxes in the first 
and third panels, respectively. We summarize the stacking re- 
sults in Table [1] As was the case in the radio, the millimeter 
and submillimeter fluxes seem to correlate with the 4.5 fim 
flux. Unfortunately, the S/N here does not allow us to mean- 
ingfully break down the sample to further study this property. 
When compared to the entire 4.5 /im sample (open downward 
triangles in FigurefTOli. the KIERO selection picks up brighter 
objects, on average, which is also similar to the case in the 
radio. 

We stacked the millimeter and submillimeter fluxes in 
the photometric redshift subsample. At both 850 /im and 
1100 /im, the mean fluxes from z > 3 sources are higher 
than those from z < 3 sources. It is well know that the ob- 
served flux in the millimeter and submillimeter is not a strong 
function of redshift at z ~ 1 to 10 for a fixed luminosity 
( lBlain&Longairi[T993h . The weighted-combined 850 and 
1100 /im result indicates that KIEROs at z > 3 are 1.5 ±0.4 
times more luminous than KIEROs at z < 3 are, on average. 
This is consistent with the result from the radio stacking. 
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As in the radio stacking analyses, here we also investigate 
the contribution from faint sources. We removed 1100 fim 
data above 3 mJy and 850 /im data above 6 mJy (dotted lines 
in Figure [Toll, and repeated all the stacking procedures. The 
results are listed in the lower half of Table [T] The general 
trend here is very similar to that in the radio, but with worse 
S/N. 

We compared the millimeter and submillimeter and 1.4 
GHz stacked fluxes and obtained rough redshift estimates. 
This "millimetric redshift" method was first published 
by ICarilli & Y un (1999) and independently developed by 
iBarger. Cowie. & Richards (2000). Here we adopt ed the for- 
mula derived by iBarger. Cowie. & Richards! (120001) . 



:0.98(5850/5l.4GHz)' 



0.26 



1. 



(1) 



In deriving this equation, Barger et al. assumed a radio spec- 
tral index of a = 0.8 and a submillimeter dust emissivity of 
(3 = 1.0, and they used Arp 220 for flux normalization. Us- 
ing our mean 850 fim flux, our mean 1100 /im flux extrap- 
olated to 850 fim assuming /3 = 1.0, and our mean 1 .4 GHz 
flux for z < 3 .5, we find an error-weighted millimetric redshift 
of z ~ 2.4. Using the stacked fluxes on undetected sources, 
we find a consistent redshift of 2.6. Despite the large un- 
certainty of this method, the results are remarkably close to 
what would be expected from the photometric redshift distri- 
bution in Figure |5j a), which is unlikely to be a coincidence. 
The agreement between these two entirely independent red- 
shift estimates strongly suggests that our photometric redshift 
estimates and stacked fluxes are fairly reliable. 

5.3. X-Ray Properties 

S even of the 196 KIERO s are in the 2 Ms Chandra cata- 
log ([Alexander et al.ll2003h . These include the two spectro- 
scopic KIEROs at z = 1.790 and 2.578. If we assume z > 2 
for t he remaining five and an X-ra y photon index of F = 1.8 
(e.g. lBarger. Cowie. & Wang|2007l) . then all seven have either 
rest-frame hard or soft X-ray luminosities exceeding 10"^^ erg 
s"' and hence are X-ray AGNs. One of the seven sources is 
the very bright 627 /iJy radio source (source 58 in Figure |7]i. 
This extremely large radio flux and its Fnoo/jm upper limit 
(1 mJy) cannot be explained by any known starburst SED at 
z > 0.5. Thus, its radio emission is almost certainly powered 
entirely by an AGN. This justifies the exclusion of this source 
in our radio stacking analysis, since our main goal is to mea- 
sure SFRs. 

Unlike the radio image, the 2 Ms Chandra observations do 
not probe into the starburst regime at z > 2, making it difficult 
to address whether the remaining 189 KIEROs contain signif- 
icant X-ray AGN activity. At z = 2 and z = 3, a 10^^ erg s"' 
X-ray luminosity corresponds to X-ray fluxes of 4.7 x 10"'^ 
and 1.7 X 10"'^ erg s"' cm"^, respectively, for F = 1.8. The 
sensitivity of the 2 Ms observations do not reach such flux 
limits, especially in the hard X-ray band (e.g., see error bars 
in FiguresfTTIandfTSIi. Thus, KIEROs undetected by Chandra 
may stifl host > 10'*^ erg s"' X-ray AGNs. 

We carried out stacking analyses on the ha rd (2-8 keV) and 
soft (0.5-2.0 keV) unfiltered X-ray images of lAlexander et al.l 
(|20C)3|) to understand the average X-ray properties of KIEROs. 
We measured X-ray fluxes using aperture diameters approxi- 
mately 2 X the PSF FWHM. Since the FWHM of the Chandra 
images vary with distance from the field center, we adjusted 
the size of the aperture accordingly. We measured the back- 
ground locally after masking the detected objects. We com- 
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Fig. 11.— Hard X-ray (2-8 keV) flux vs. 4.5 fj.m flux for KIEROs. The 
bottom panel shows the X-ray fluxes of individual sources. Sources enclosed 
in large squai'es are AGNs identified in Section |6] The top panel shows the 
stacked X-ray fluxes in each 0.5 dex 4.5 fim flux bin, only for sources be- 
low the dotted line, which roughly corresponds to 5(t. Solid squai'es are the 
mean fluxes for KIEROs. The exact values are also shown near the symbols. 
Triangles show the mean fluxes of the entire 4.5 /im sample in the same flux 
ranges. In the two lower 4.5 ^m flux bins, both the KIERO sample and the 
4.5 /^m selected sample are not detected. 
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Fig. 12.— Same as Figure [T7] but in the soft X-ray (0.5-2.0 keV). The 
middle panel is a blowup of the bottom one. 

pared our X-ray flux measurements with those given in the 
2 Ms catalog and found excellent agreement. We estimated 
the errors (caused by confusion and shot noise) in the stacked 
fluxes with Monte Carlo simulations of random sources, as 
we did in the longer wavebands. 

The stacked hard and soft X-ray fluxes for all the 196 
KIEROs are {Fhx) = 2.9±0.9 x lO"'^ and (Fsx) = 8.1 ± 1.4 x 
10""* erg s"' cm"^. The hard X-ray value is close to the 
flux for X-ray AGNs, and the soft X-ray value falls signifi- 
cantly below. Again, to understand if these stacked fluxes are 
dominated by small numbers of bright X-ray sources, we ap- 
plied flux cutoffs that are approximately 5 a in both hard and 
soft bands, indicated by dotted lines in Figures [TT| and [121 
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We imposed the same cutoff in the Monte Carlo simulations, 
which allows us to greatly reduce errors caused by bright X- 
ray sources that get close to our targets by chance. Further- 
more, we stacked the X-ray fluxes of the KlEROs according 
to their 4.5 fim fluxes, as we did in the longer wavebands. The 
results are presented in the top panels of FiguresfTTIandfTSl 

In Figure [TT] the stacked hard X-ray fluxes are consis- 
tent with AGNs at > 3 /iJy. This suggests a signif- 
icant AGN contribution in the E( 5 brightest members of 
KlEROs. On the other hand, KlEROs with < 3 fiJy 
show no detectable hard X-ray fluxes, and all the stacked soft 
X-ray fluxes in Figure[T2]are below the flux required for X-ray 
AGNs. This suggest that the majority of KlEROs do not host 
X-ray AGNs. In the following sections, we will further iden- 
tify AGNs in the radio and mid-lR (MIR) and exclude these 
AGNs in our analyses of star formation properties. 

5.4. MIPS 24 Properties 

The Spitzer GOODS Legacy Program DRIh- data release 
provided a 24 pm catalog obt ained with the Spit zer Multiband 
Imaging Photometer (MIPS: iRieke et al.ll2004l) . The catalog 
has a 24 pm flux limit of 80 pJy and has 1041 sources in our 
IRAC region. Among the 1041 sources, 18 are KlEROs. The 
fraction of KlEROs in the DRIh- 24 pm sample is thus 1.7%, 
lower than the fraction in the radio (5.1% of the 1 .4 GHz sam- 
ple of iMorrison et al.ll20Toi Section lU). However, the 80 
flux limit used in the DRIh- catalog is shallower than the 
true depth of the image. To search for fainter MIPS detected 
KlEROs, we first made a dee p source extraction in the MIPS 
image with SExtractor (iBertin. & Arnoutslll996h . We looked 
for matches between KlEROs and faint 24 pm sources within 
2" radii (approximately 1/3 of the PSF FWHM). We then in- 
spected all the matches by eye to remove unreliable matches 
caused by nearby bright sources and noise spikes. In partic- 
ular, there are KlEROs residing in crowded regions with sev- 
eral blended MIPS sources. We removed such KlEROs unless 
their are clearly associated with local 24 pm peaks. The above 
procedure increases the number of MIPS detected KlEROs to 
27. Simple aperture photometry indicates that this pushes the 
24 fim detection limit to slightly less than 50 fiJy, correspond- 
ing to approximately 5 cr. After this, the fraction of 24 pm 
detections among the sample of 196 KlEROs is 14%, compa- 
rable to the fraction of radio detected KlEROs, which is 11%. 
This similarity suggests that radio and 24 i^im observations are 
roughly equally effective in picking up KlEROs. 

There is, however, an additional complication in the 24 pm 
case. At z > 1, multiple broad PAH emission and silicate ab- 
sorption features enter the MIPS 24 i^im band, producing both 
positive and negative selection biases that are strong func- 
tions of redshift. These strongly affect the observed 24 pm 
fluxes of KlEROs, since they are all at z > 1. Because of 
this, using 24 pm fluxes to quantify the st ar formation activity 
in KlEROs is highly un certain (see, e.g.. lMurphy et al.ll2009l : 
iRodighiero et alJl2010h and strongly depends on the assumed 
SED models. On the other hand, the SEDs of galaxies in the 
radio are simple power laws with very similar spectral indices, 
and the radio-FlR correlation is fairly insensitive to dust tem- 
perature. We thus only relied on the radio data to study the 
SFRs of KlEROs. 

We did not attempt to probe deeper at 24 pm with stack- 
ing analyses given the large uncertainties in interpreting the 
results. Instead, we stacked the radio, millimeter, and submil- 
limeter fluxes of the MIPS detected and undetected KlEROs, 
to see whether the stacked fluxes are dominated by the few 



24 pm detected sources and hence whether there are also 
passive objects in the KIERO selection. First, the 27 MIPS 
detected KlEROs have mean 1.4 GHz, 1100 ^m, and 850 
pm fluxes of 42.1 ± 1.7 pJy (after excluding the > 600 pJy 
bright radio source), 1.94 ±0.22 mJy, and 6. 15 ±0.68 mJy, 
respectively. All these values are approximately 4x higher 
than the mean fluxes of all KlEROs listed in Table [T] This 
is not surprising since we expect the 24 pm detections to 
be very luminous at high redshift. The remaining non-MlPS 
KlEROs have mean 1.4 GHz, 1 100 ^m, and 850 pm fluxes of 
5.0±0.7 ^Jy, 0.27±0.09mJy, and0.45±0.31 mJy, respec- 
tively. The results are still significant, except at 850 pm where 
the sample size is small (72). Furthermore, if we exclude the 
^^1.4 GHz > 20 pJy sources from our radio stacking analyses, 
as we did in Section 15.11 then the mean radio flux of the re- 
maining 157 non-MlPS sources is 2.3 ±0.5 piy. This is com- 
parable to the result of 2.6 ± 0.5 piy for all F14 ghz < 20 piy 
sources. We conclude that even sources undetected at 24 pm 
have significant radio, millimeter, and submillimeter emis- 
sion. We do not find evidence for a passive population based 
on the 24 pm to 1.4 GHz data. 

6. AGN FRACTION 

In Section 15.31 we fo und seven KlEROs de tected in the 
Chandra 2 Ms catalog ([Alexander et al.l l2003h . All seven 
sources are X-ray AGNs with soft or hard X-ray luminosities 
greater than 10"*" erg s"'. To obtain a more complete iden- 
tification of X-ray AGNs, we measured the X-ray fluxes of 
KlEROs as we did in the X-ray stacking analyses. We found 
four additional sources with > 3(7 X-ray fluxes and with X- 
ray luminosities exceeding lO'*^ erg s"'. Therefore, there are 
1 1 X-ray AGNs in our KIERO sample^ 

In the MIR, our pre vious work jBarger. Cowie. & Wand 
l2008l: IWang et alj2010h had sh own that thelRACAGN selec- 
tion techniques of lLacv et al] ( 120041) and lStern et al. I ( l2005h 
will either be incomplete or highly contaminated by normal 
galaxies in a deep sample like t his. Here we adopt the MIR 
power-law selection technique (iDonley et alJl2007l) . We re- 
quire that the objects have to be detected in all four IRAC 
bands, and we performed power-law fitting to their IRAC 
SEDs. Among the 196 KlEROs, 54 are bright enough to be 
detected in all IRAC bands, and 13 have SEDs that can be 
well fitted ix' probability > 90%) by red power laws. All of 
them have red IRAC spectral slopes with a > 0.7, consistent 
with them being AGNs. One of the 13 is detected by Chandra 
and is also an X-ray AGN. 

In the radio, we estimated the rest-frame 1.4 GHz power 
of radio detected KlEROs using their photometric redshifts 
and assuming a synchrotron spectral slope of a = 0.8. We 
found two KlEROs with unusually large 1 .4 GHz radio pow- 
ers of > 10^^ W Hz"'. Neither is a MIR power-law AGN. 
One so urce is #58 in Figure |7] an X-ray AGN discussed in 
Section 15.31 Its photometric redshift is very likely to be a 
catastrophic failure caused by its unusual SED, but this does 
not affect its AGN identification. The other source is #70 
in Figure |2l Its photometric redshift fitting appears reason- 
able, and thus the estimated rest-frame radio power is reliable. 
Its large radio power (2.1 x 10^^ W Hz"') requires a SFR of 
> 10"* Mq yr"' (see next section), which is not supported by 
its infrared properties (e.g., it is undetected at 24 pm). We 
thus conclude that its radio emission is powered by an AGN. 

The above X-ray, MIR, and radio selection identify 23 
AGNs, which is 12% of the KlEROs in the GOODS-N. This 
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fraction is a lower limit, since AGNs with lower X-ray, MIR, 
or radio luminosities would not be selected. On the other 
hand, all 23 AGNs are detected in the IRAC bands. Thus, the 
AGN fraction of the 54 members of the IRAC-bright subsam- 
ple is ^ 40%. This sets an upper limit on the AGN fraction of 
the entire KIERO sam ple since we expect to have more AGNs 
in luminous sources. iMessias et al.l (120101) studied AGNs in 
high-redshift extremely red galaxy populations using radio. 
X-ray, and MIR criteria. They found an AGN fraction of up to 
^ 30%, among which the majority are obscured MIR power- 
law AGNs. Their results are consistent with ours, given the 
small sample size here. 

To see if the above AGN identification picks up X-ray 
AGNs that are below the Chandra detection limit, we repeat 
the hard X-ray stacking analyses in Figure [TT] without all the 
above AGNs. In the 4.5 /im flux ranges of > 10 /iJy and 3-10 
/iJy, the stacked hard X-ray fluxes decrease substantially to 
4.0 ±2.5 X lO"''^ and 0.3 ±2.1 x lO"'^ erg s"' cm'^, respec- 
tively. Both values are not statistically significant. Therefore 
the above X-ray, MIR, and radio AGN identification is fairly 
effective in removing potential X-ray AGNs. 

7. STAR FORMATION RATE 

The SFRs of the KIEROs can be inferred from their total 
IR luminosities. Among all the flux measurements in Ta- 
ble [T] the radio fluxes provide the most robust estimates of 
the FIR luminosities. This is because the radio spectral slopes 
of normal galaxies are very similar and there is a tight cor- 
relation between radio and IR luminosity of norm al galax- 
ies, re gardless of dust temperature (see the review of lCondonl 
119921) . The latest FIR observations suggest tha t this corre- 
lation also holds for high-redshift galaxies (e.g.. lIvison et alj 
|2010a b; Bourne et al. 201 1). An additional reason for using 
the radio is the higher S/N in the stacked radio fluxes. 

A possible uncertainty of using radio fluxes for SFR esti- 
mates is the contribution from AGNs. In the previous sec- 
tion, we identified AGNs using X-ray, MIR, and radio data. 
This allows us to exclude AGNs from our sample. Further- 
more, most of the KIEROs have radio fluxes well below 100 
/iJy. At z ~ 3 this flux corresponds to a radio power of 
6 X 10""* W Hz"'. In the local universe, AGNs dominate the 
radio luminosity funct ion at Pi 4 ghz > 10^^ W Hz"' (e.g., 
iMauch & Sadler! I2007h . On the other hand, at z > 1, star- 
bursts start to dominate sources around lO^"* W Hz"', and the 
luminosity function of AGNs falls substantially below that of 
starbu rsting galaxies at lower radio powers (e.g.. lCowie et alj 
120041) . Thus, it is plausible that the contributions from AGNs 
to the radio fluxes of our KIERO sample are negligible. Below 
we will show this is indeed the case using AGNs identified in 
the previous section. 

We now derive the conversion between radio flux and SFR. 
The radio power can be calculated from 

PlA GHz = 4TTdlil+zr-'FiA GHz, (2) 

where di is the luminosity distance and a is the radio spectral 
slope. Here we assume a = 0.8. To calibrate the conversion 
between radio power an d SFR, we assume the radio-FlR cor- 
relation (ICondonll 19921) 

^ " U.75XOT ) -''^ [whF^ ) 

a nd the can onical q value of 2.3. Using the SED templates 
m ISflvaetal..(.1998.) , we found that for a broad range of FIR 



SEDs the total IR luminosity conventionally defined between 
8 and 1000 /im is approximately 2x that of the 40-120 fim 
FIR luminosity. With this conversion factor and th e standard 
conversion between total IR luminosity and SFR in lKennicutl 
(fT998h . 

SFR (Mq yr"') = 1 .7 x 10"'°L,R(io), (4) 

we derived 

SFR(M0yr"') = 7.83x 10"^(l+zr"' ( -^V ( . 

\MpcJ V MJy / 

, ^ 

This agrees with the conversion in lYun. Reddv. & Condor 
f2001) to within 10% and is ^ 20% higher than that inlBell 

In Table |2] we present the star formation properties of the 
KIEROs based on the stacked radio fluxes. The upper half of 
the table is for all KIEROs, and the lower half is for KIEROs 
without X-ray, MIR, and radio AGNs. It can be seen that ex- 
cluding AGNs only slightly decreases the mean radio fluxes, 
except in the highest redshift bin, where the stacked radio 
fluxes are significantly biased by the two radio AGNs. This 
supports the argument that most of the KIEROs have radio 
emission powered by star formation. Below we will limit our 
discussion to the non-AGN results. 

For the sources with redshifts, the SFRs and IR luminosi- 
ties are si milar to those o f local ultraluminous IR galaxies 
(ULIRGs; ISanders & MFabeL 1996) such as Arp 220. They 
are also close to the lower end of typical dusty galaxies found 
by millimeter and submillimeter single-dish surveys (Ljr > 
IO'^Lq). This is what we expected for the KIERO color se- 
lection. It is difficult to estimate the SFRs of the KIEROs 
without redshifts. If their redshift distribution is similar to 
that of KIEROs with redshifts, then they would be roughly 
an order of magnitude less luminous, or luminous IR galaxies 
(LIRGs). On the other hand, it is possible that these Ks and 
IRAC faint sources are at higher redshifts. A mean radio flux 
of 1.9 fiJy from z ^ 3-4 would imply a SFR of ~ 10^ Mq yr"' 
and a ULIRG luminosity. We can thus safely conclude that 
most of the KIEROs are at least LIRGs at high redshifts, and 
many of them are starbursting ULIRGs. 

We investigated whether there is a correlation between stel- 
lar mass and SFR. We did this only in the z = 2-4 range be- 
cause of the sample size. For z = 2-3, we divided the non- 
AGN KIEROs into three groups according to their stellar 
masses: < lO'^M© (7 sources), 1O"'-1O"M0 (11 sources), 
and > 10"Mq (7 sources). Their stacked radio fluxes are, re- 
spectively, 12.7 ±3.3, 5.6 ±2.6, and 23.3 ±3.3 /iJy. Forz= 3- 
4, we divided them into two groups: < 10"Mq (6 sources), 
and > 1O"M0 (9 sources). Their stacked radio fluxes are, 
respectively, 19.2 ±2.9 and 18.4 ±3.5 /iJy. There is not a 
convincing trend here, except for fluctuations likely caused 
by the small sample sizes. Therefore, even the most massive 
KIEROs in our sample are still actively forming stars, and 
there is no evidence for a significant number of massive and 
passive galaxies in the KIERO population. This is quite dif- 
ferent than other high-redshift red galaxies. We believe this is 
because our KIERO color selection avoids the 4000 A Balmer 
break (a signature of old stars) and primarily picks up galax- 
ies with continuous red spectral slopes (a signature of dust 
extinction). 

We estimated the contributions of KIEROs to the cosmic 
star formation, i.e., their SFRD. Normally, one would like to 
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TABLE 2 

Star Formation Properties of KIEROs 



Redshift 


^galaxy 


{FlA GHz) 




(SFR) 


SFRD 






(/iJy) 


(IQi^Lo) 


(Mo yr-') 


(Mq yr-l Mpc-3) 


All Sources 


z~ 1.5-2 


10 


5l.0±2.7 


3.4 


580 


0.016 


z~2-3 


27 


13.4±1.7 


2.2 


370 


0.013 


z~3-4 


23 


18.5±1.8 


6.7 


1150 


0.036 


z~ 4-5.5 


10 


31.3±2.7 


20.8 


3532 


0.053 


unknown 


120 


2.5 ±0.8 








Non-AGNs 


z~ 1.5-2 


9 


46.5 ±2.9 


3.1 


530 


0.013 


z~2-3 


22 


14.5 ± 1.9 


2.4 


400 


0.011 


z~3-4 


15 


18.9±2.2 


6.9 


1170 


0.024 


z~ 4-5.5 


7 


7.3 ±3.3 


4.8 


820 


0.009 


unknown 


115 


1.9 ±0.8 









apply the 1 /Vmax method (ISchmidtlll968h for calculating vol- 
ume. Here, since most of our sources are not individually 
detected in the radio, applying such a method is effectively 
introducing additional weighting according to their Ks or 4.5 
l^im properties in the radio stacking analyses. It is unclear to 
us whether this would bias the results. Therefore, we simply 
adopt total comoving volumes in the four redshift bins in Ta- 
ble |2] for our survey area (0.06 deg^). The derived SFRD for 
the sources with redshifts are listed in the table and shown 
in Figure [13] (solid squares). These values are lower limits 
for the following reasons: (1) they are not corrected for the 
substantial incompleteness in the Ks and 4.5 fim bands; (2) 
they do not include star formation in the KIEROs identified 
as AGNs; and (3) they do not include contributions from the 
KIEROs without redshifts. If we assume that the redshift dis- 
tribution of the KIEROs without redshifts is similar to that of 
the KIEROs with redshifts, then the SFRD values in Tabled 
would increase by 20% in all four redshift bins. On the other 
hand, if we assume that the KIEROs without redshifts are, on 
average, from higher redshifts, then their contribution to the 
SFRD would be substantial. For example, if we assume that 
they are at z > 3, then the SFRDs at z > 3 would increase by 
60% (open squares in Figure [iJTl. 

The SFRD values we found ar e comparable to thos e for 
bright SMGs at similar redshifts ( IChapman et alj|2005h . but 
we note that not all the bright SMGs are KIEROs. In Wang 
et al. (2006), we found a nearly flat SFRD of - 0.06-0.15 
Mq yr"' Mpc"-' at z = 1-4, from our radio and 850 fim stack- 
ing analyses of a large 1.6 and 3.6 fim selected NIR sam- 
ple (dashed line in Figure [13] while the dotted line shows 
the same SFRD maximally corrected for incompleteness in 
the z = 1-3 range). Naively speaking, the KIEROs should 
be a subsample of the larger NIR sample of Wang et al. 
(2006). However, the CLEAN-like procedure employed in 
WIO for this work provides many more 3.6 and 4.5 fim de- 
tected sources and better photometry. By comparing the list of 
KIEROs with our 3.6 fim catalog compiled in 2006, we found 
that 45% of the KIEROs in the SCUBA area were missed by 
Wang et al. (2006). The missed KIEROs have a broad range 
of 3.6 /im fluxes and the same IRAC image is used in both 
works, so this is not an effect of different detection limits. 
Because of this, comparing the SFRD here with that in Wang 
et al. (2006) is not straightforward. 

We also compare our SFRDs with recent extinc- 
tion corrected r est-fra me UV determinations from LBGs 
(iBouwens et al j [ 20091 open diamonds in Figure [T3l 
21' 




iReddv & Steidelll200^ 



open 
open triangles). 



At 



in 

z > 3, 



the SFRD 



3 4 

Redshift 

Fig. 13. — Comoving SFRD. Solid squares aie SFRDs derived from our 
radio stacking analysis on KIEROs with redshift information. Open squares 
further include the contributions from KIEROs without redshifts, by assum- 
ing z > 3 for these Ks and IRAC faint KIEROs (see text for details). Open 
diamonds and triangles are extinction conected rest-frame UV results from 
[Bpuwens et al. (2009) and Reddv & Steidel (2009), respectively, for LBGs 
integrated to 0.04 L*. The dashed line shows the mean SFRD in Wang et 
al. (2006) at z = 1^, derived using a radio stacking analysis on a large NIR 
selected sample. It is approximately flat. The dotted line shows the SFRD in 
Wang et al. (2006) at z = 1-3 assuming the maximum possible completeness 
conection in the submillimeter. 

of the KIEROs is similar in magnitude to the SFRD of the 
extinction corrected LBGs. At z < 3, the KIERO SFRD con- 
tribution becomes only ^ 10% of that of LBGs. This might 
be explained by the fact that KIEROs do not include all SMGs 
and therefore they only account for part of the dusty star for- 
mation. However, even the radio/submillimeter SFRD de- 
termined by stacking on a large NIR sample by Wang et al. 
(2006; dashed line in Figure [T3j is significantly lower than 
the SFRD of the extinction corrected LBGs. 

The extinction corrected SFRD of LBGs seems to be un- 
usually large at z > 2, compared to submillimeter and KIERO 
SFRDs. In Table [1] we showed that optically faint KIEROs 
are brighter in the radio, millimeter, and submillimeter than 
KIEROs detected by ACS. In Section [O] we find that 
KIEROs have very little overlap with LBGs. One important 
finding in Wang et al. (2006) is that most of the submillimeter 
EBL arises from galaxies with intermediate rest-frame colors, 
i.e., galaxies not as blue as LBGs. The combination of these 
points suggests that a substantial amount of high-redshift star 
formation cannot be traced by rest-frame UV emission, and 
that the total SFRD should be the sum of the SFRD from 
LBGs and the SFRD from dusty sources. This would give 
a very large total SFRD at z > 2, which raises the question of 
whether the large extinction corrected LBG SFRD is overes- 
timated. 
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ICarilli et all ( 120081) stacked radio fluxes of z - 3 LBGs in 
the COSMOS field and found that the radio inferred mean 
SFR of these LBGs is only 1.8 times the extinction uncor- 
rected mean UV SFR. They pointed out that this is much 
smaller than the factor of ~ 5 gen erally adopted for high- 
redshift LBGs. Furthermore. iHo et al . (2010) could not detect 
the radio signal from z ^ 4 LBGs with an extremely deep ra- 
dio stacking analysis in the GOODS-N and S fields. Their 
inferred radio SFR of z ^ 4 LBGs i s only consiste n t with 
the extinction corrected UV SFR in iBouwens et al.l (120091) 
if a 2 (T radio flux upper limit is adopted. This result also 
suggests a discrepancy between the radio SFR and the dust- 
corrected UV SFR. The reason for this discre pancy is un- 
clear. In addition to the possibilities discussed in lCarilli et alj 
( 120081) . another possibili ty is th e correla tion between FIR flux 
and UV spectral slope (iMeure r. Heck man. & Calzettil Il999h 
assumed in many high-redshift LBG studies. Re cent observa- 
tions s how substantial sca t ter in this correlation (ICortese et al 
l2006t iHowell et all I20T0I: iWiiesinghe et al.1 iMlT and this 
may contribute to the uncertainty in the UV extinction correc- 
tion. This issue will require further investigation, especially 
once direct measurements of the FIR luminosities of LBGs 
with sensitive instruments such as ALMA and Herschel be- 
come available. 

8. KIERO AND OTHER HIGH-REDSHIFT POPULATIONS 
8.1. IRAC Selected EROs 

lYan et al.l (l2004l) used the ACS z' (F850LP) image and the 
IRAC 3.6 ^m image in the HUDF and found 17 IRAC se- 
lected ext remely red objects (lEROs) with z'-3.6 /^jm > 3.25 
(also see rYanll2008h . Among our 196 KIEROs, 123 have 
z' detections at > 3cr, and only 27 of them are above the 
lERO color cut. If we adopt the 3(t limits for z' undetected 
KIEROs, then the number of lEROs increases to 68. Thus, 
> 1/3 of KIEROs are also lEROs. Among the 17 lEROs in 
Yan et al., five satisfy the KIERO color cut, including one 
Ks undetected lERO. Thus, approximately 1/3 of lEROs are 
KIEROs. The cumulative density of lEROs with Ks < 25 is 
~ 1.4 arcmin"^, more than 2x higher than that of KIEROs 
even after a Ks based completeness correction. From these 
points of view, KIEROs and lEROs are two different popula- 
tions with roughly 30% overlap. 

On the other hand, the KIEROs are remarkably similar to 
lEROs in many other properties. First, the 17 lEROs in Yan 
et al. have redshifts between 1.6 and 3.6, 1 1 of them between 
2 and 3. This is almost identical to the redshift distribu- 
tion for KIEROs (Figure |5ji). Secondly, the SEDs of lEROs 
mostly consist of strong NIR components from massive pop- 
ulations of old stars, plus secondary rest-frame UV compo- 
nents from ongoing star formation. This is also very similar 
to the SEDs of KIEROs (Figure|2]and SectionlH, except that 
some KIEROs are too faint in the optical for the rest-frame 
UV components to be detected. The stellar masses derived 
from the SED fitting for lEROs and KIEROs are also simi- 
lar, approximately 10^^-10"^Mq. In Section|2] we showed 
that 14 out of the 28 identified SMGs in the GOODS-N are 
KIEROs. If we include /-band upper limits, then 22 out of 
the 28 SMGs are lEROs. This fraction is even higher than 
KIEROs. ^ 

Furthermore, iMessias et alj ( 120101) performed radio stack- 
ing analyses on several red galaxy populations including 
lEROs. They found a mean radio flux of 8.2 ±1.9 from 96 
lEROs at z = 2-3 after excluding 3 sources detected at > 3a, 



where 1 cr is ^ 14-17 fiJy in their observations. With a similar 
^ 50 /iJy cut in our stacking of z = 2-3 KIEROs, the stacked 
radio flux is 7.5 ± 1 .4 /iJy, quite comparable to lEROs at sim- 
ilar redshifts. The only difference here is that nearly 10% of 
KIEROs (2 out of 27) are above the 50 /iJy flux cut, but only 
3% of lEROs in Messias et a l. (201Q) are above this cut, in- 
dicating that the KIERO population contains more ULIRGs. 
However, within the numerical uncertainties, these two per- 
centages could be consistent. 

Why two such similar populations do not overlap much in 
terms of color properties and number densities is an interest- 
ing question. In the color-color diagram in Figure [14] we 
see that the lERO selection and the KIERO selection over- 
lap at high redshift on dusty sources. The difference is in 
the low-redshift end. In the lower-right part of the diagram, 
early-type galaxies (including those not forming stars) enter 
the lERO selection at z > 1 .5 or so. Such galaxies do not en- 
ter the KIERO selection at z < 3. The upper-left part of the 
diagram is more interesting, i.e., galaxies that are KIEROs 
but not lEROs. Such galaxies have blue z'-3.6 i-im colors 
that are consistent with starbursting galaxies at z > 6, simply 
based on the color tracks. These K IERO s are not z > 6 galax- 
ies. Instead, as described in Section l431 they have two distinct 
stellar populations. Their massive and reddened stellar popu- 
lations make them redder than the KIERO criterion, and the 
unobscured parts of their ongoing starbursts make them bluer 
than the lERO criterion. In short, lEROs contain more early- 
type galaxies, and KIEROs contain more dusty sources with 
ongoing star formation. We believe this is a consequence of 
the fact that the KIERO selection avoids the Balmer break but 
targets the red spectral slopes caused by strong dust extinc- 
tion. 




2.5log(F3.6/F,-) 

Fig. 14. — z', Ks, and IRAC color-color diagram. The vertical and horizon- 
tal lines indicate the color criteria for lEROs and KIEROs, respectively. Color 
curves show the tracks of elliptical galaxies (thin red cu rves, Cole man et al. 
1980) a nd starburst galaxies ( thick blue curves, Kinney e t al.lI996l) . reddened 
with the lCalzetti et al.l 120001) extinction law with Av=Q (solid), 1 .0 (dotted), 
2.0 (dashed), and 3.0 (dash-dotted). Integer redshifts are indicated by crosses 
along the tracks, starting from i = in the lower-left and ending at z = 6. 
Solid diamonds show KIEROs detected at both Ks and z' ■ Open diamonds 
with an'ows show 1 a z' — 3.6 /im limits of z' -undetected KIEROs. Dots in 
the background are field galaxies that are detected at z, Ks, and 3.6 ^m at 

> 5 CT. 
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8.2. Distant Red Galaxies 

In additi on to lEROs, dist ant red galaxies (DRGs, J - 
Ks > 1 .35: iFranx et al.l 120031) are another interesting pop- 
ulation with which to compare. DRGs have redshifts 
roughly between 2 and 4 (e.g.. Ivan Dokkum et ail llOOl 
iFoster Schreiber et al.1 l2004t iReddv et al.l l2005h . similar to 
KIEROs and lEROs. We therefore expect a substantial over- 
lap between the DRG and KIERO populations. Unfortunately, 
there is not yet a /-band image in the GOODS-N that is suf- 
ficiently deep and wide to verify this. The CFHT WlRCam 
y-band image that we mentioned in Section |4] covers the en- 
tire GOODS-N with an rms depth of 0.09 fi]y (cf. 0.12 /xJy 
at Ks). Of the 104 Ks detected KIEROs, only 37 are detected 
in this /-band image, and of those, only six have J -Ks col- 
ors redder than 1.35. An additional 15 7 undetected KIEROs 
have 7-band flux upper limits (3<t) satisfying the DRG color. 
Thus, at least 52 of the 104 Ks detected KIEROs are DRGs, 
and the DRG fraction among the Ks detected KIEROs can be 
anywhere between 50% and 70%. 

There exists another /-band image around the HDF-N taken 
with the 8.2-m Subaru Telescope (Kajisa wa et al. 200 6). It is 
substantially deeper than the CFHT ima ge. We carried ou t 
an independent reduction of this image (iWang etal.ll2007h . 
which covers ~ 26 arcmin^ with an rms depth of ~ 0.03 /iJy. 
There are 16 Ks detected KIEROs in this image. Five of them 
are detected in this deep / image but all of them do not satisfy 
the DRG criterion. Five of the remaining six 7-band unde- 
tected sources have J upper limits that satisfy the DRG crite- 
rion. Thus, the fraction of Ks detected KIEROs that are DRGs 
is anywhere between 30% and 70% based on this narrower 
and deeper / image. 

Finally, in the extremely deep HST NICMOS Fl lOW image 
of the HDF-N (e.g., Dickinson et al..2000; rms - 0.15 /iJy), 
there are only three KIEROs. Only one of them is detected in 
the /Tj-band and it is not a DRG. However, this image is too 
narrow for us to draw a meaningful conclusion. 

In terms of number density, the cumulative density of DRGs 
at Ks < 24 is > 2 arcmin"^ dLabbe et al.i ,2003: Grazian et al.l 
l2006bl:lKaiisawa et al.l2006h . significantly higher than that for 
KIEROs (0.9 arcmin^ at our Ks and 4.5 iim limits). However, 
there is tentative evidenc e that the number den sity of DRGs 
starts to drop at Ks > 24 (iKajisawa et al.| |2006l), while that of 
KIEROs still seems to be increasing at this magnitude (Sec- 
tionO. If both observed trends are real, KIEROs are not just 
the faint-end tail of the DRGs. Based on the above results, we 
conclude that most DRGs are not KIEROs. On the other hand, 
a substantial fraction (perhaps > 50%) of KIEROs are DRGs, 
but not all. In Figure[T5]we show the J, Ks, and 4.5 /^tm color- 
color diagram. It shows that the difference between DRGs 
and KIEROs is somewhat similar to that between lEROs and 
KIEROs. 

In the SMG population, the fraction of DRGs is high. 
Among the 28 identified SMGs in the GOODS-N, 15 are 
DRGs. This fraction is comp arable to the KIERO fraction in 
SMGs. iKnudsen et al.l (120051) studied the SFRs of DRGs in a 
lensing cluster field with a submillimeter stacking analysis at 
850 ^m. For DRGs with K < 24.4 they found a mean 850 fim 
flux of 0.93 mJy and a mean SFR of 127 Mq yr"' (both cor- 
rected for lensing amplification). Based on this, the typical 
SFRs of KIEROs are > 2x higher than those of DRGs, de- 
spite b eing fainter in the NIR. On the radio side. lMessias et al.l 
(120101) performed radio stacking analyses on DRGs. They 
found a mean radio flux of 6.1 ± 1.6 /iJy from 152 DRGs at 




-10 12 3 4 
2.5log(FK3/Fj) 

Fig. 15. — /, Ks, and 4.5 fim color-color diagram. The vertical and hor- 
izontal lines indicate the color criteri a for DRGs and KIEROs. respectively. 
Curves are the same as those in Figure[T4] Solid diamonds show KIEROs de- 
tected in both Ks and /-band WIRCam images. Open diamonds with arrows 
show 1 cr J —Ks limits of 7-undetected KIEROs. Dots in the background are 
field galaxies that are detected at /, Ks, and 4.5 fim at > 5 cr. 

z = 2-3 after excluding 3 sources detected at 3ct. As men- 
tioned in Section 18.11 a similar stacking of z = 2-3 KIEROs 
gives a comparable 7.5 /iJy radio flux, but the KIERO popu- 
lation has a higher ULIRG fraction. 

8.3. Lyman Break Galaxies 

To find LBGs in the 112 members of the ACS -detected 
subsample of KIERO s, we adopted the selection criteria of 
iBeckwith et"aI1 (l2006h for galaxies in the H UDF. The selec- 
tion cr iteria are also similar to those used by iBouwens et al.l 
(120071) . We found 15 ^-dropouts, 4 v-dropouts, and 5 i- 
dropouts. This shows that most of the KIEROs are not LBGs, 
even for the subsample that is detected in the optical. The rea- 
son for this is that the spectral slopes are too red to be selected 
by the dropout criteria. Furthermo re, the number of LB Gs 
(- 1700 /7-dropouts in GOODS-N; 'Bo uwens et all 120071) is 
much higher than the number of KIEROs (197 in GOODS - 
N). All these results imply that KIEROs and LBGs are two 
distinct populations. 

It is important to realize that KIEROs are the most inten- 
sive starbursting galaxies at z > 2, second only to the bright- 
est SMGs. Their high SFRs (Section |7]i imply strong rest- 
frame UV radiation. Without absorption, their redshifted op- 
tical fluxes would be roughly 1.6-2 /iJy at 0.5-1 /im. How- 
ever, even in the ACS-detected KIERO subsample, most of 
them have optical fluxes lower than this (e.g.. Fig. |7]i. The 
factor of almost 10 in the rest-frame UV dust attenuation in 
KIEROs is much larger than the generally a ccepted extinc- 
tion corrections of 3 x to 5 x in LBGs (e.g., IBouwens et al.l 
l2007t IReddv etaHIlM IBouwens et al.ll2009l) . 

8.4. BzK Galaxies 

The BzK selection (iDaddi et al.ll2004l) is an effective tech- 
nique for finding 14 < z < 2.5 galaxies. Since this redshift 
range overlaps with the lower-end of the redshift distribution 
of KIEROs, it is interesting to see where KIEROs fall on 
the BzK diagram (Figure [TSIl. A difficulty here is that most 
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KIEROs are optically faint, and we can only put a small sub- 
sample of KIEROs on the BzK diagram. Only 21 KIEROs are 
detected in all three of the b, z', and Ks bands (squares), 14 of 
which are starburst BzK galaxies (sBzK). Sixteen KIEROs are 
detected in both the z' and Ks bands (triangles), 1 1 of which 
have b-z' lower limits consistent with sBzK. The non-BzK 
KIEROs either have photometric redshifts that are too high or 
too low for the BzK selection (#15, 35, 53, and 55 in Figure|7]i 
or have unusually blue rest-frame UV and optical SEDs (#44 
and 63 in Figure|7]i that barely miss the BzK selection. 

We conclude that roughly 2/3 of the optically bright 
KIEROs are sBzK galaxies. However, the requirement of b 
and z'-band data makes the BzK selection less effective for 
extremely optically faint populations like the KIEROs. Fi- 
nally, among the 28 identified SMGs, 13 are sBzK galaxies. 
This fraction is si milar to the KIERO fr actio n in SMGs. How- 
ever, a s shown by lPannella et al.l ( l2009h and lKurczynski et"an 
(1201 Ih . the BzK selection picks up much more "normal" star 
forming galaxies with SFR 10-100 Mq/ji. The SMG frac- 
tion in BzK galaxies is much lower than the SMG fraction in 
KIEROs. 
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Fig. 16. — b, z' , and Kg color-color diagram. The black diagonal and 
horizontal lines are the selection crit eria for sBzK and pBzK galaxies. Curves 
are the same as those in Figure [14] except that we added Scd-type galaxies 
(green curves, Coleman et al. 1980). Solid diamonds are KIEROs detected 
at Ks, z', and b bands. Open diamonds with aiTows show 1 a z' — b limits 
of fo-undetected KIEROs. Dots in the background are field galaxies that are 
detected in all three of the b, z' , and Ks bands at > 5 cr. Squares are KIEROs 
that were detected in all three bands at > 3 cr. Triangles show the lower limits 
for the KIEROs that are undetected in the fo-band. 



9. DISCUSSION 

9.1. Contributions to the Extragalactic Background Light 

Our original goal of selecting extremely red sources from 
the Ks and IRAC 4.6 /im bands was to see if we could se- 
lect high-redshift dusty sources in the NIR. In this paper 
(Section |7| we showed that KIEROs are mostly LIRGs and 
ULIRGs at z > 2. We found that KIEROs are optically faint. 
Their red optical SEDs imply that they will not be picked out 
as LBGs, even if they are detected in the optical. All of the re- 
sults indicate that we have successfully selected a population 
of galaxies that are too dusty to be selected by existing optical 



and NIR surveys and too faint to be systematically detected 
with existing MIR, millimeter, and submillimeter surveys. 

An interesting question is whether the KIERO sample in- 
cludes the majority of ULIRGs at z > 2. At 850 /im and 
1100 /im, KIEROs have an integrated surface brightness of 
4.47 ±0.88 and 1.64 ±0.26 Jy deg"- (Table [B, respectively. 
The EBL measured by the COBE satellite has strengths of 3 1- 
44 Jy deg-2 at 850 /im and 18- 25 Jy deg" - at 1 100 /im. These 
were measured by two groups dPuget et a l. 1996; Fixsen et"an 
119981) . and the range reflects the uncertainties in removing 
foreground contamination. Thus, KIEROs contribute approx- 
imately 10% to the EBL at these wavelengths and do not seem 
to represent the majority of sources that give rise to the mil- 
limeter and submillimeter EBL. 

On the other hand, there is evidence that the majority 
of the millimeter and submillimeter EBL (es pecially at > 
850 um) arises at z < 2 (Wang et al. 2006; ISerieant et alj 
120081: iMarsden et alJl2009l) . and thus that the KIERO selec- 
tion may still pick up a large fraction of z > 2 dusty sources. 
For example, Wang et al. (2006) found an 850 /im EBL con- 
tribution of ^ 16 Jy deg"- from a sample of 1.6 and 3.6 /im 
selected sources at z < 2. This leaves at most 15-28 Jy deg"^ 
to sources at z > 2, of which ^ 15%-27% can be attributed to 
the KIERO sample here. The actual fraction should be even 
larger considering the fact that the KIERO selection is still 
severely limited by completeness at Ks and 4.5 /jm. To deter- 
mine whether or not KIEROs can fully account for the EBL 
that arises at z > 2, we will need both a better determination 
of the KIERO number counts and a better EBL measurement. 

9.2. Dust Hidden Star Formation 

In Section 18.31 we showed that the dust attenuation of the 
rest-frame UV radiation of KIEROs is roughly a factor of 10 
or even higher It is now known that some of the most lu- 
minous millimeter or submillimeter selected galaxies can be 
entir ely hidden by dust in th e opti cal and even in th e NIR 
(e.g.. lWang. Barger. & Co ^¥ 2009'; ICowie et al.l l2009'). The 
work presented here shows that there are still extremely extin- 
guished galaxies when we go roughly an order of magnitude 
fainter in the total IR luminosity. Our KIERO color selection 
picks up such galaxies at z ^ 1 .5-5, but the Ks and IRAC im- 
ages become less sensitive to galaxies at z > 3. With a narrow 
selection window between z 2 and 4 and the completeness 
limits in the Ks and the 4.5 /im bands, we already pick up 
roughly 10% of the total EBL at 850 /im and 1 100 /im. The 
total fraction of background that arises from extremely extin- 
guished galaxies from z = to > 6 should be even higher than 
this. We thus expect a significant fraction (^ 10%) of cosmic 
star formation to be entirely hidden from deep optical obser- 
vations. 

9.3. Nature of KIEROs 

The KIERO color selection is meant to pick up objects with 
red spectral slopes over broad wavelength ranges and to avoid 
the 4000 A Balmer break. This makes it sensitive to dusty 
starbursting galaxies. The fact that the KIEROs are selected 
near the peak of their stellar SEDs (i.e., the rest-frame 1.6 /im 
bump) also means that it is sensitive to massive systems. In 
SectionlHwe showed that most KIEROs have stellar masses of 
10'" to 10'- Mq. In Section|7]we showed that most KIEROs 
are ULIRGs with Lir > lO'^L© and are actively forming stars. 
The masses of KIEROs are very simi lar to those of bright 
850 /im selected SMGs at z > 2 (e.g., iDve et al...2008.) . On 
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the other hand, the mean 850 flux of KIEROs (1.44 mJy) 
is below the confusion limit of single-disk telescopes. Its cor- 
responding IR luminosity (assuming a dust SED similar to 
Arp 220) are a few times lower than those of most SMGs se- 
lected at 850 fim by single-dish telescopes. This suggests that 
KIEROs are massive galaxies undergoing slightly less inten- 
sive starbursts compared to bright SMGs. 

In Section|2]we found that half of the identified SMGs in the 
GOODS-N are KIEROs. The fraction of SMGs that are DRGs 
is comparable (Sections 18.2b . and the fraction of SMGs that 
are lEROs is significantly higher (^ 80%, Section lOl ). How- 
ever, this does not mean that all of these color selections are 
equally efficient in picking up dusty sources. The lERO and 
DRG selections pick u p substantial numbe rs of passive and 
z<2 galaxies (see also. lMessias et al.]|2010l) . However, in our 
radio stacking analyses of KIEROs, we did not see evidence 
for a significant subpopulation of passive KIEROs. This is 
true even after we introdu ce a 24 /im selection to our radio 
stacking analyses (Section 15.4b . It is especially remarkable 
that approximately 1/3 of KIEROs are not detected by ACS. 
Such UV-faint sources would be considered as "passive" by 
many studies, but these KIEROs are even brighter in the ra- 
dio than UV-bright KIEROs, indicating that they are more ac- 
tive in star formation. This distinguishes the KIERO selection 
from the lERO and DRG selections. We believe the KIERO 
selection is more effective in picking up dusty starburst galax- 
ies and contains fewer passive galaxies. 

In this paper we provide estimates of stellar masses and 
SFRs of KIEROs. Another crucial parameter for understand- 
ing the evolution is gas mass, which is not yet established 
for this new population of galaxies. Given the great simi- 
larity (SFRs, dust obscuration, and stellar masses) between 
KIEROs and bright SMGs, it is likely that the two populations 
are similarly massive in molecular gas mass. Recent observa- 
tions of high-y CO transit ions of SMGs found molecular gas 
masses of a few 10 ' °M,t. (iNeri et al.l l200l iGreve et aljr2005l: 
iTacconi et all 120061 l2008h . If KIEROs are similarly rich in 
molecular gas, they can sustain their ultraluminous starburst 
phase for '^100 Myr. In other words, they can significantly 
increase their stellar masses in less than 10% of their Hubble 
time. 

To summarize, the KIERO selection seems to pick up a 
sample of massive and dusty z > 2 galaxies. They are form- 
ing stars actively. This is consistent with massive galaxies at 
z> 2 that are still rapidly building up their mass. The prop- 
erties of KIEROs are similar to their more luminous coun- 
terparts (SMGs), as well as to high-redshift massive galaxies 
selected in the IRAC bands (e.g., lEROs and DRGs; see also 
iMancini et al.ll20d9l) . In each of these various selections we 
are likely seeing different regions of the luminosity function 
or different stages in the formation of the most massive galax- 
ies at z > 2. Eventually these galaxies will become quiescent 
and passively evolve into the "downsizing" era. However, at 
high redshifts, most of them seem to be still active. Because 
of the dust obscuration, common for the most intensive star- 
forming systems, a complete view of such massive and active 
galaxies requires observations at both long and short wave- 
lengths. 

10. SUMMARY 

In order to find high-redshift, faint, dusty sources that may 
give rise to a significant fraction of the submillimeter EBL, 
we selected 196 KIEROs with extremely red colors of Ks- 
4.5 /im> 1 .6 in the 0.06 deg^ GOODS-N region. The selected 



KIEROs have a range of 4.5 /im fluxes, mostly between 1 and 
10 ^l]y. Of the 196 KIEROs, 104 KIEROs are detected in 
the Ks band, and the slope of the Ks number counts for these 
sources is steeper than that of Ks selected galaxies at similar 
flux levels, likely because KIEROs are at redshifts higher than 
Ks selected galaxies. The counts also still seem to be rising at 
the Ks detection limit of - 0.3 /iJy. Roughly 2/3s of KIEROs 
are detected in the optical by HST ACS. However, only 1/5 
of the ACS-detected KIEROs are LBGs. The remaining ACS 
detected KIEROs all have optical continua that are too red to 
pass the LBG selection. 

We performed a photometric redshift analysis on a Ks- 
bright subsample of 76 KIEROs. We found a redshift distri- 
bution between z ^ 1 .5-5, with roughly 70% at z 2-A. This 
redshift distribution is quite similar to that of bright SMGs, 
but with a bias to high redshift. Our photometric redshift 
analysis also shows large stellar masses of 1O'"-1O'^M0. We 
found that a significant fraction of KIEROs have SEDs that 
cannot be well fitted by simple stellar populations. They have 
massive old stellar components that are luminous in the NIR 
and ongoing starbursts that are detected in the optical. 

We found that roughly half of the identified millimeter 
and submillimeter sources in the GOODS-N are selected as 
KIEROs, but such sources only represent 7% of the selected 
KIEROs. To study the radio, millimeter, and submillimeter 
properties of the remaining KIEROs, we carried out stacking 
analyses. We found that, on average, KIEROs are brighter 
than normal 4.5 fim selected galaxies by factors of ^ 2 and 
^ 6-10 in the radio and in the millimeter and submillime- 
ter, respectively. We also found that KIEROs undetected by 
ACS are ~ 3-4 times brighter in all these wavebands than 
KIEROs detected by ACS. We repeated the stacking analyses 
on KIEROs undetected in the radio, millimeter, and submil- 
limeter. These faint KIEROs still possess detectable stack- 
ing signal, showing that their properties are not dominated by 
small numbers of bright sources. We also did not find evi- 
dence of passive KIEROs based on the 24 /im to radio data. 

We identified AGNs in the KIERO sample with large X-ray 
and radio luminosities, and with MIR SEDs that are feature- 
less red power laws. We identified 19 such sources, of which 
2/3 are obscured AGNs identified in the MIR. We also found 
an AGN fraction of up to 1/3. 

Using the local radio-FIR correlation of normal galaxies 
and the stacked radio fluxes of non-AGN KIEROs, we found 
that they are LIRGs and ULIRGs. The results are consistent 
with the KIEROs being a population of z > 2 ultraluminous 
dusty galaxies that are slightly fainter than SMGs in the mil- 
limeter and submillimeter and undetected by single-dish sur- 
veys. Their large IR luminosities imply high SFRs of ~ 300- 
1200 Mq yr"\ and therefore a population of massive z > 2 
galaxies in rapid formation. We showed that KIEROs con- 
tribute roughly 10% to the total EBL at 850 and 1 100 /im and 
up to 30% of the EBL at these wavelengths which arises from 
galaxies at z > 2. Given that KIEROs and LBGs have very 
little overlap and that optically faint KIEROs are brighter in 
the radio, millimeter and submillimeter, these results indicate 
that a significant fraction 10%) of cosmic star formation is 
hidden by dust and cannot be traced by rest-frame UV emis- 
sion. 
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